ACCELERATED  LIFE  TESTING  AND  RELIABILITY  PREDICTION 

by 

SAM  RUSTOMJEE  DARUVALLA 

B,  Sc,  (Physics-Mathematics),  University  of  Bombay 
Bombay,  India,  1958 


13/ 


A  MASTER'S  THESIS 


submitted  in  partial  fulfillment  of  the 


requirements  for  the  degree 


MASTER  OF  SCIENCE 


Department  of  Electrical  Engineering 


KANSAS  STATE  UNIVERSITY 
Manhattan,  Kansas 


1965 


Approved  by 
Major  Profes'sor 


TV  ^^ 

PP77  TABLE  OF  CONTENTS 

'^^       -I- 

Page 

INTRODUCTION  1 

DEFINITIONS  4 

LIFE-TESTING  TKHNIQUES 6 

AXELERATED  LIFE-TESTING  TECHNIQUES ■. 11 

Step-Stress  Aging  12 

Estimation  of  Acceleration  Factors 15 

Problem  of  Modeling 16 

RELIABILITY  MODELS  OF  A  COMPLEX  SYSTEM 18 

Exponential  Model  , 18 

Weibull  Model 25 

Markovian  Model 33 

Worst-Case  Design  Model 36 

CORRELATION  TB3HNIQUES 40 

RELIABILITY  PREDICTION  FOR  A  SIMPLE  AMPLIFIER  51 

CONCLUSIONS 60 

APPENDIX  I  62 

ACKNOWLEDGEMENTS 63 

BIBLIOGRAPHY  64 


INTRODUCTION 

During  the  last  two  decades,  the  trend  towards  automation  and 
precision  products  in  industry  and  military  defense  systems  has  led 
to  ever-increasing  complexities  in  equipment  design  and  the  consequent 
difficulty  in  the  prediction  of  equipment  reliability.  The  over-all 
reliability  of  a  complex  system  with  no  redundant  parts  is  the  product  of 
the  reliability  of  its  components,  as  given  by  the  formula 


n 

R     ,,  =  n  R.  (1) 

over-all  •—,     ^ 
1-1 


where  R.  is  the  reliability  of  the  i-th  component.  The  equation  indicates 
that  each  individual  component  reduces  the  over-all  reliability  of  such  a 
system  by  its  ovm  reliability  factor.  A  simple  system  consisting  of,  say, 
100  components,  each  having  a  99  per  cent  reliability,  will  have  an  over- 
all reliability  of  only  36.5  per  cent.  A  more  complex  system  consisting  of 
400  components,  ivith  99  per  cent  individual  component  reliability,  would 
have  an  over-all  reliability  of  only  3  per  cent. 

Hov;ever,  a  modern  commercial  airliner,  with  all  its  complexities,  is 
known  to  be  more  reliable  than  the  above  discussion  would  suggest.  This  is 
so  because  only  a  few  of  its  components,  most  of  the.;',  structural,  are  really 
vital  in  the  sense  that  failure  of  any  one  of  them  v.'ill  cause  a  total  loss 
of  the  aircraft.  Many  other  components,  particularly  the  electronic  com- 
ponents, are  not  vital  since  even  in  the  event  of  their  failure,  the  air- 
craft can  still  be  operated  safely.  On  the  other  hand,  each  component  of 
the  guided  missile  system  is  vital  since  failure  of  any  one  of  them  will 


cause  the  complete  missile  system  to  fail.  Hence,  the  components  used  in 
a  missile  must  be  made  perhaps  four  orders  of  magnitude  (or  ten  thousand 
times)  more  reliable  than  commercial  components.  In  other  v;ords,  the 
achievement  of  near  absolute  reliability  of  components  used  in  missile 
would  be  a  desirable  goal  to  be  specified  in  their  design  and  manufacture. 
In  order  to  specify  that  a  component  has  a  reliability  of  99.99  per  cent 
per  1000  hours,  with  99  per  cent  confidence,  it  v;ould  be  necessary  to  run 
a  life  test  of  a  representative  sample  for  more  than  46,000,000  unit-hours 
until  the  first  failure  occurred. 

It  is  clearly  seen  that  in  the  case  of  life  testing  under  simulated 
normal  operating  conditions,  a  decrease  in  the  duration  of  the  tests  is 
accompanied  by  a  corresponding  increase  in  the  number  of  components  tested. 
However,  it  is  possible  to  cause  a  small  sample  to  fail  in  a  short  length  of 
time  by  increasing  the  environmental  stress  level  to  many  times  that  under 
normal  operating  conditions.  Thus,  the  failure-rate  distribution  curves 
under  accelerated  life-testing  conditions  are  obtained  vjith   considerable 
saving  in  time  and  cost. 

The  question  then  arises  v/iiether  it  is  possible  to  predict  the  relia- 
bility of  the  component  under  normal  operating  conditions  from  the  results 
of  accelerated  life-tests.  In  those  cases  where  it  is  possible  to  establish 
a  statistical  correlation  between  the  failure-rate  distributions  under  the 
two  different  sets  of  environmental  conditions,  the  component  reliability 
under  normal  operating  conditions  can  be  predicted  from  the  accelerated  life- 
test  data.  The  results  of  accelerated  life-tests  must  be  applied  with  great 
care  since  the  accelerated  conditions  are  far  removed  from  the  noimal  operat- 
ing conditions.  Despite  the  drawbacks  encountered,  accelerated  life-testing 
is  becoming  increasingly  important  since  it  alone  offers  a  solution  to  the 


dilemma  of  testing  highly  reliable  components  used  in  military  defense 
systems. 

A  complex  missile  system  costs  millions  of  dollars,  and  hence,  the 
simulated  and  accelerated  life-testing  of  complete  missiles  is  not  mone- 
tarily feasible.  Besides,  a  missile,  once  fired,  cannot  be  retrieved 
and  hence,  in  the  event  that  a  missile  fails  to  hit  its  target,  the  causes 
of  failure  cannot  be  analyzed  and  corrective  procedures  applied.  This  then 
sets  up  a  new  problem  which  is  the  prediction  of  reliability  of  a  complex 
system  from  the  data  on  accelerated  life-testing  of  its  components. 

An  analytical  treatment  of  various  reliability  models  for  a  complex 
system,  namely,  the  exponential  model,  the  Weibull  model,  the  Markovian 
model  and  the  worst-case  design  model,  is  presented.  It  is  shov;n  that  each 
of  these  models  can  explain  certain  aspects  of  system  reliability,  but  be- 
cause of  the  various  assumptions  involved,  none  of  these  can  full  explain 
the  anomalies  sometimes  observed  when  the  reliability  predictions  are 
subjected  to  operational  verification.  The  reliability  of  a  simple  ampli- 
fier circuit  is  predicted  using  the  exponential  model  which  assumes  the 
system  failure-rate  is  the  sum  of  the  individual  component  failure-rates, 
provided  the  latter  remain  constant  over  the  useful  life  of  the  system. 

Finally,  the  technique  of  correlating  the  accelerated  failure-rate 
data  to  the  failure-rate  under  normal  operating  conditions  is  discussed, 
wherein  it  is  assumed  that  the  failure  modes  of  the  system  are  the  same 
for  various  stress  level  combinations  of  different  environments.  A  more 
general  correlation  technique  for  correlating  the  data  in  case  of  a  non- 
exponential  behavior  of  failure-rate  with  environment  is  also  presented. 
The  latter  technique  can  be  applied  even  in  cases  where  the  condition  of 
same  failure  modes  under  the  two  sets  of  environmental  conditions  may  not 


hold,  provided  the  acceleration  curve  begins  to  take  up  its  general  shape 
within  the  region  of  accelerated  environmental  stress  levels, 

DEFINITIONS 

Reliability.  It  is  the  probability  of  a  device  performing  its  purpose 
adequately  for  the  period  of  time  intended  under  the  operating  conditions 
encountered.  Adequate  performance  implies  that  a  prespecified  criteria  is 
satisfied  by  the  device  under  test.  The  period  of  time  and  the  operating 
conditions  affect  the  system  performance  and  they  are  included  as  part  of 
reliability  specifications. 

System.  A  combination  of  subsystems  which  is  capable  of  operation  by 
itself;  for  example,  a  guided  missile. 

Subsystem.  A  group  of  components,  combined  and  packaged  in  one  housing 
to  perform  a  particular  set  of  operations;  for  example,  a  radio  transmitter, 
a  servo  amplifier,  etc. 

Component.  An  item,  such  as  a  resistor,  capacitor,  electron  tube,  gyro, 
etc.  which  is  not  normally  subject  to  further  disassembly. 

Mean-Time-to-Failure.  The  ratio  of  the  total  test  time  of  a  system  or 
a  component  to  the  total  number  of  failures.  For  a  life-test  experiment  in 
which  n  components  are  put  on  test  and  the  test  is  terminated  after  the  first 
r  failures,  the  mean-tLme-to-failure  of  the  components  is 


it.   -^  (n  -  r)t^ 
9  =  -ill (2) 


where  t^  is  the  time  counted  from  the  beginning  of  the  test  to  the  i-th  failure. 
Equation  (2)  assumes  that  the  components  that  fail  are  not  replaced  by  fresh 


components.  In  case  of  testing  v;ith  replacement  the  mean-time-to- failure 
of  the  components  is  given  by 


nt 

e  =  — ^  (3) 


If  the  testing  is  discontinued  after  a  certain  fixed  time  T,  regardless  of 
the  number  of  components  that  failed,  then 


e  =  -l~  (4) 

T 


where  r_  is  the  total  number  of  components  failing  in  time  T. 

Failure-Rate.  '  It  is  the  reciprocal  of  mean-time-to- failure.  It  is  de- 
noted by  A  and  is  often  expressed  as  failures  for  1000  hours  or  as  per  cent 
failure-rate  per  1000  hours, 

^  =  -f  (5) 

Failure-Modes.  It  is  the  physical  state  of  the  component  at  the 
instant  of  failure.  For  example,  two  of  the  possible  failure  modes  of 
a  resistor  are  "open"  and  "short"  states.  For  a  capacitor,  dielectric 
breakdovm  is  one  of  the  failure  modes. 

Catastrophic  Failures.  These  are  failures  which  cause  a  normally 
operating  system  to  suddenly  become  completely  inoperative,  for  example, 
a  random  "open"  occurring  in  a  wire  resistor  after  several  hundred  hours 
of  operation.  This  type  of  failure  is  usually  caused  by  chance  and  there- 


fore  cannot  be  predicted  in  advance  for  a  particular  time  and  a  specific 
part,  without  the  knowledge  of  the  failure-rate  distribution. 

Degradation  (Drift)  Failures.  These  are  the  type  of  failures  which 
occur  gradually  due  to  the  change  of  a  parameter  viith  time,  for  example 
a  change  in  the  ohmic  value  of  a  resistor  may  cause  excessive  hum  or 
frequency  shift  in  a  radio  receiver. 

Wear-Out  Failures.  These  are  failures  v;hich  can  be  predicted  on  the 
basis  of  a  knovm  v/earout  characteristic,  for  example,  the  v;ear  of  the 
brushes  of  an  electric  motor. 

Secondary  Failures.  A  secondary  failure  occurs  as  a  result  of  a 
primary  failure.  If  a  resistor  in  an  electronic  circuit  is  shorted,  caus- 
ing an  excessive  drain  on  a  tube;  the  tube  failure  is  considered  as  secondary 
failure  because  it  was  due  to  the  primary  resistor  failure. 

LIFE-TESTING  TECHNIQUES 

In  recent  years,  considerable  interest  has  been  stimulated  about  the 
analysis  of  life-test  data,  A  life-test  is  conducted  by  subjecting  one  or 
more  identical  components  to  a  given  set  of  operating  conditions  and  noting 
the  number  of  hours  of  satisfactory  performance,  that  is,  the  time-to-failure 
of  each  component.  Failure  is  defined  as  the  state  in  which  the  component 
no  longer  performs  satisfactorily.  The  techniques  discussed  by  Zelen  (1959) 
involve  the  analysis  of  failure-rate  data  for  a  range  of  environmental  con- 
ditions. Oftentimes,  life-testing  procedures  involve  the  acceleration  of 
the  failure  of  components  by  testing  at  over-stress  environmental  conditions. 
The  accelerated  life-test  data  are  then  extrapolated  to  normal  operating 
stress  levels.  The  various  criteria,  which  must  be  satisfied  if  such  extra- 
polation are  to  be  valid  and  meaningful,  have  been  outlined. 


It  is  assumed  that  the  time-to-failure  follows  an  exponential  distri- 
bution having  the  probability  density  function,  given  belov/ 


1  exp  [-  I  (t-Y)  ]         t  s  Y 
:(t)  =  {  Q       '  (6) 

^  0  t  <  Y 


where  G  is  the  mean-time-to-failure  or  mean  life  of  the  components.  The 
reciprocal  of  mean  life  is  the  mean  failure-rate.  The  above  distribution  is 
called  a  tv;o-parameter  exponential  distribution  in  (9,  y)»  where  the  parameter 
Y  has  the  property  that  no  failures  occur  before  time  y  expires. 

The  discussion  is  limited  to  two  different  environments  E-,  and  E^  with 
p  and  q  stress  levels,  respectively.  There  are  pq  different  combinations 
of  environments  and  for  each  combination,  n  components  are  simultaneously 
placed  on  test.  The  experiments  are  terminated  when  exactly  r  of  the  com- 
ponents on  test  have  failed.  The  r  failure  times,  for  each  treatment  combi- 
nation, are  denoted  by  t,  ^  t„  ^  ...  ^  t  . 

A  particular  treatment  combination  is  denoted  by  the  symbol  ij  Vi/here  i 
and  j  refer  to  the  stress  levels  of  E  and  E^,  respectively.  The  effects 
of  various  treatment  combinations  on  the  mean-time-to-failures,  8.  .  of  the 
components,  irrespective  of  the  value  of  y- •  ^re  analyzed.  Let  G. .  be  de- 
fined as 

e..  =  mp.q.k..     '  ■'^^;"'   ;      _  (7) 


where  p.,  q.  are  positive  constants  representing  the  contributions  from  the 
stress  level  i  of  environment  E  and  the  level  j  of  E^,  respectively.  The 
quantity  k. .  is  a  constant  which  depends  on  both  environments,  and  m  is  a 


constant  common  to  all  treatment  combinations. 

If,  for  all  treatment  combinations,  p.q.k^.  =  1  then,  9^.  =  m  and 
6.  .  is  independent  of  the  stress  levels  of  the  environments  of  E  and  E^. 

If  p.k. .  =  1  (or  q.k.  .  =  l)  then,  9^.  =  mq.  (or  9^.  =  mp^)  and  9^. 

depends  on  the  stress  levels  of  E^  alone  (or  stress  levels  of  E  alone). 

If  k.  .  =  1,  G. .  =  mp.q.  and  hence,  the  effects  of  environment  E,  can  be 

evaluated  independently  of  the  effects  of  environments  E„  and  vice-versa. 

Thus,  v;hen  9.  .  =  mp.q.,  it  is  seen  that  the  effects  of  the  tv;o 
'       ij    ^x^3 ' 

environments  are  completely  multiplicative.  With  such  a  model,  accelerated 
life-testing  is  feasible  because  the  stress  level  of  any  one  environment 
can  be  chosen  unusually  high  in  order  to  increase  the  component  failure- 
rate,  keeping  the  stress  level  of  the  other  environment  at  the  normally 
operating  value. 

Consider  an  experiment  conducted  over  all  possible  treatment  combi- 
nations of  the  stress  levels  of  the  environments  E  and  E„,  and  let 

t,  £:  to  ^  ...  ^  t  be  the  time-to-failure  of  the  first  r  components  out  of 
i    ^         r 

the  initial  n  components  put  on  test.  The  miniraum  variance  unbiased  esti- 
mate of  G. .  is 


e   =  r(t  -  ti)  -i-  (n  -  r)  (tr  -  ti)  ,   . 

ij  (r  -  1)  ^  ^ 


where 

t  =  4-   E  t  (9) 

i=l 


The  maximum  likelihood  estimators  of  the  various  parameters 
m,  p. J  q.  and  k. .  are 


1/, 


A 

/    P          C       A       ) 

m  = 

TT         TT      0-  . 

i=ij=i    ') 

pq 


^ 


A 

P.- 


1 

"TT 

m 


rr  e.  . 


,Vc 


i  =  1,2,3  ... 


A 

k. 


ij 


A 

'  A'  A TT- 

m  p.  Q. 


(10) 


A 
m 


TT  e.  _ 
1=1  "■' 


'^ 


3   =  1,2,3  ... 


Tnus,  if  it  is  desired  to  predict  the  value  of  9. .  for 
say;  it  is  required  to  evaluate 


=  3  and  j  =  2, 


A  A   A  A 
ra  P3  q^   k32 


This  should  agree  with  the  observed  value  of  G„p,  except  for  rounding  errors. 


10 


In  some  applications,  where  the  parameter  y  is  taken  to  be  equal  to 
zero,  the  minimum  variance  unbiased  estimate  of  9.  .  is  modified  to 


A   _  t  -(n-r)t^  ^^) 


The  above  techniques  for  predicting  9. .  were  based  on  the  assumption 
that  the  failure  times  have  an  exponential  distribution.  This  leads  to 
the  question  as  to  v/hat  procedures  should  be  used  if  the  underlying  distri- 
bution of  time-to-failure  is  other  than  exponential.  Fortunately,  all  results 
for  the  exponential  distribution  can  also  be  used  for  distributions  having 
the  probability  density  function  of  the  form 


,  h(t)  exp  (  -f  h(x)  dx)  t  ^  o 

'(*'  =  {  0         °  t  <  0         (12) 


The  function  h(t)  is  a  non-negative  function  of  the  time-to-failure  t.  The 
exponential  and  V/eibull  distributions  are  special  cases  of  the  general  family 
of  distributions  described  by  (12).  If  then,  H(t)  follows  an  exponential 
distribution  with  mean-time-to-failure  equal  to  9.  It  is  always  possible  to 
transform  the  time-to-failure  to  a  quantity  H(t)  which  will  follow  an  expo- 
nential distribution.  Hence,  all  formulae  based  upon  the  exponential  as- 
sumption can  be  used  for  the  transformed  variate  H(t). 

For  example,  in  the  case  of  random  failures,  the  probability  of  failure 
is  independent  of  the  time  that  the  component  has  been  performing  satis- 
factorily. In  such  a  case,  the  quantity  h(t),  the  probability  that  the 
component  will  fail  in  the  next  interval  t  and  t  +dt  is  independent  of  t. 


11 


that  is 


h(t)=-^  (14) 


and  hence 


f(t)  =  -i-  exp  (  -  i  )  (15) 


the  exponential  probability  density  function. 

For  the  case  of  wearout.  failures,  h(t)  is  an  increasing  function  of 
time,  say 


h(t)  =  -3£1.  ■  (16) 

6 


This  yields  the  Weibull  probability  density  function  for  v;hich 

H(t)  =  t^  (17) 

ACCELERATED  LIFE-TESTING  TECHNIQUES 

It  is  shovm,  in  Appendix  I,  that  in  order  to  predict  with  99  per  cent 
confidence  that  a  component  under  test  has  a  reliability  of  99.99  per  cent 
per  1000  hours,  it  is  necessary  to  run  a  life-test  of  a  statistically  repre- 
sentative sample  for  more  than  40,000,000  unit-hours  until  the  first  failure 
occurs.  The  present  day  life-testing  techniques  although  theoretically 
feasible  are  not  economical,  and  therefore,  accelerated  life-testing  is  re- 
sorted to.  This  section  covers  accelerated  life-testing  techniques  and  the 


reliability  prediction  at  normal  operating  stresses  from  accelerated  life- 
test  data. 

Accelerated  life-testing  is  a  misnomer  because  it  is  the  component 
failure-rate,  and  not  its  life,  which  is  accelerated.  The  environmental 
conditions,  in  such  tests,  are  more  severe  than  those  under  normal  conditions. 
Such  conditions  are  achieved  by  increasing  the  ambient  temperature,  the  ap- 
plied voltage,  etc.  or  combinations  of  these  factors.  This  invariably  re- 
duces the  duration  of  the  life  test  by  causing  early  failures  of  the  components 
on  test,  • .'  " 

The  following  three  methods  represent  the  commonly  used  techniques  in 
the  field  of  accelerated  life-testing:  '  ■ 

1.  Accelerated  life-testing  by  step-stress  aging. 

2.  Accelerated  life-testing  by  estimating  acceleration  factors. 

3.  Accelerated  life-testing  as  a  problem  of  modeling. 

Each  method  yields  significant  information  and  each  has  its  limitations. 

Step-stress  aging.  An  accelerated  life-testing  program  is  used  to 
obtain  failure-rate  acceleration  curves,  as  shov/n  in  Figure  1.  For  a  particu- 
lar environmental  test  condition,  say  S  ,  a  series  of  experiments  are  conducted 
to  obtain  component  failure-rate  data,  using  a  fresh  sample  from  the  same  sta- 
tistical population  each  time.  The  failure-rate  distribution  curve  for  this 
particular  environmental  stress  condition  is  obtained,  as  shown.  The  experi- 
ments are  then  repeated  for  different  environmental  stress  levels,  each  more 
severe  than  the  normal  operating  environmental  stress  level.  The  mean  of  the 
failure-rate  distribution  is  a  function  of  the  environmental  stress  level  and 
an  acceleration  curve  is  dravm  through  these  points.  If  the  failure-rate 
distributions  have  similar  distributions  at  each  stress  level,  it  indicates 
that  the  failure  modes  for  the  components  are  the  same  for  the  different 


13 


w 
u 

8 

o 
o 
o 


u 
a> 

a. 

0) 
•4-> 
(0 
^1 
I 


(0 
Mh 

+> 

C 
u 


Acceleration 
curve 


Failure-rate 
distribution 


normal 


S. 

1 


Stress  level 


Figure  1:  Accelerated  failure-rate  curve  for  various  stress  levels 
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Figure  2:  Accelerated  stress  curve  for  various  failure-rates 


environmental  stress  levels.  In  such  a  case,  true  acceleration  has  been 
achieved  and  it  is  possible  to  predict  the  failure-rate  at  normally  operating 
stress  levels  by  extrapolation  of  the  accelerating  curve. 

Another  method  of  obtaining  the  acceleration  curve  is  to  fix  the 
failure-rate  of  the  component  and  determine  the  amount  of  stress  necessary 
to  cause  component  failure,  as  shown  in  Figure  2.  The  above  techniques  of 
step-stress  aging  are  particularly  suitable  for  use  in  early  stages  of  com- 
ponent development  to  compare  small  batches  of  components  manufactured  by 
different  methods. 

Estimating  Acceleration  Factors.  This  method  requires  the  knowledge 
of  analytical  relationship  between  the  number  of  components  and  their  time- 
to-failures.  For  example,  in  case  of  the  capacitors,  the  ratio  of  the  mean 
lives  of  two  samples  of  capacitors  is  inversely  proportional  to  the  p-th 
power  of  the  applied  voltages,  all  other  environmental  factors  remaining  the 
same;  that  is. 


=  M  '      '      ■  (IB) 


VI 


where  9^  and  82  are  mean  lives  of  the  two  batches  of  capacitors  when  the 
applied  voltages  across  their  terminals  are  v.  and  v^,  respectively.  The 
exponent  p  is  called  the  acceleration  factor  for  the  applied  voltage. 

Accelerated  testing  is  conducted  on  fresh  batches  of  capacitors  be- 
longing to  the  same  statistical  population  at  different  applied  voltages. 
The  mean  life  for  a  particular  applied  voltage  can  be  estimated  as 
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A    t,  +  t^  +  ...  +  t  +(n-r)t  (19) 

fl  -   1 4 r      _    T_ 


A 

where  9   is  the  estimated  mean  life 


t.   is  the  time-to-failure  of  the  i-th  component 


r   is  the  number  of  failures  after  which  the  experiment 
is  terminated,  and       -  r  • 

n   is  the  number  of  components  initially  on  test. 
The  negative  slope  of  the  log-log  plot  of  mean  life  versus  applied  voltage 
is  the  estimate  of  p.  It  is  then  possible  to  evaluate  the  mean  life  of 
the  capacitor  at  rated  applied  voltage,  with  the  use  of  equation  (l8). 

The  accuracy  of  this  technique  depends  on  the  nianber  of  failures  and 
on  the  assumption  that  the  modes  of  component  failure  remain  the  same  under 
different  stress  levels  of  the  combination  of  environments. 

A  Problem  of  Modeling.  It  is  assumed  that  the  principles  in- 
volved in  the  prediction  of  performance  characteristics  of  a  component  from 
that  of  its  model  can  also  be  applied  to  the  prediction  of  component  reli- 
ability from  the  tests  associated  with  its  model.  The  basic  assumption  in 
the  use  of  accelerated  testing  is  that  the  probability  of  failure  is  a 
function  of  test  duration  and  stress.  In  this  section,  the  measure  of  time 
is  not  the  total  number  of  operating  hours,  but  instead,  the  actual  time 
spent  under  high  stress  levels  which  cause  signification  degradation  of 
the  system. 

The  probability  of  failure  of  a  component  is  assumed  to  be  dependent 
only  on  the  cumulative  effects  of  high  stress  levels  over  the  entire  test 
duration.  Suppose  that  the  ratio,  A,  of  the  time  for  which  the  component 
is  subjected  to  excessive  stresses  to  the  total  test  time,  is  known.  If 
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the  characteristic  time  T  be  defined  as  the  number  of  operating  hours 
required  for  the  component  to  accumulate  1  hour  of  overstress  degradation, 
then 

By  increasing  the  abuse  ratio  to  A'  for  the  model,  the  characteristic  time 
T^  for  the  model  can  be  decreased.  For  a  good  analysis,  the  abuse  ratio  a' 
for  the  model  should  be  unity,  that  is,  the  model  is  subjected  to  excessive 
stresses  throughout  the  test  duration.  If  T*  is  the  total  test  time  for  the 
model,  then,  the  corresponding  time  T  for  the  component  is 


(-^)  " 


(21) 


where  A   is  called  the  time-scale  factor.  Equation  (21)  is  valid  only  in 

Pi 

the  case  of  catastrophic  failures. 

Suppose  that  for  a  particular  test  conducted,  the  abuse  ratio  is  found 

to  be  0.001,  that  is,  1  hour  of  abuse  is  accumulated  for  every  1000  hours 

of  operation.  Further,  suppose  that  10  models  are  operated  simultaneously, 

each  at  an  abuse  ratio  of  0.9.  The  time-scale  factor  ^'  _  900,  so  that  1 

A 

hour  of  model  operation  is  equivalent  to  900  hours  of  component  operation. 
If,  for  the  models,  the  first  failure  occurs  after  1000  hours  of  operation, 
this  corresponds  to  1  failure  out  of  10  components,  operating  at  normal  en- 
vironmental conditions,  at  the  end  of  9  x  10^  hours. 

The  preceding  three  methods  of  accelerated  life-testing  require  that 
the  tests  should  be  relatively  simple  and  inexpensive  and  that  true  acceler- 
ation must  be  achieved.  The  last  requirement  is  somewhat  difficult  to  achieve 
as  the  accelerated  testing  conditions  are  far  removed  from  the  operating 
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conditions  resulting  in  different  modes  of  failure  for  the  component.  Hov/- 
ever,  if  the  results  are  applied  with  care,  accelerated  life-testing  tech- 
niques present  the  only  feasible  way  out  of  the  dilemma  of  reliability 
prediction  of  a  complex  system. 

RELIABILITY  MODELS  OF  A  COf^IPLEX  SYSTEM 

A  complex  system,  such  as  a  missile  control  system,  consists  of  thousands 
of  components  and  costs  millions  of  dollars  and  so  it  is  not  economically 
feasible  to  conduct  any  type  of  life  tests,  simulated  or  accelerated,  on  the 
complete  system.  In  order  to  predict  total  system  performance,  various  relia- 
bility models  are  used  which  employ  life-test  data  on  individual  components 
comprising  the  system.  In  this  section,  four  such  reliability  models,  namely, 
the  exponential  model,  the  Weibull  model,  the  Markovian  model  and  the  worst- 
case  design  model,  are  discussed  and  it  is  shov/n  that  each  of  these  models, 
while  explaining  certain  aspects  of  system  reliability,  cannot  full  explain 
the  anamolies  observed  when  reliability  predictions  are  subjected  to  oper- 
ational verification. 

The  Exponential  Model.  It  assumes  that  the  failure-rate  distri- 
bution of  any  component  can  be  approximated  by  three  different  segments,  as 
shown  in  Figure  3.  The  early  failure  region,  or  "infant  mortality"  region, 
in  which  the  component  has  a  decreasing  failure  rate,  followed  by  a  long  low 
stretch  of  chance  failure  region,  in  which  the  failure-rate  is  more  or  less 
constant,  and  then,  the  v/ear-out  region  during  which  the  failure-rate  increases. 

The  magnitude  and  length  of  the  constant  failure-rate  region,  also  called 
the  normally  operating  region,  depends  upon  the  different  environmental  stresses. 
Consider  a  component  having  a  constant  failure-rate,  X,  for  a  particular  combi- 
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Figure  3:  Failure-rate  versus  test  duration  for  the  exponential  model 
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nation  of  environmental  stresses,  being  subjected  to  life  tests.  In  a 
small  interval  of  time  (t,  t  +  At),  the  probability  of  failure  is  X   At, 
independent  of  previous  history  and  time  because  of  the  assumption  of 
constant  failure-rate.  If  R(t)  is  the  probability  that  the  component 
survives  t  time  units,  that  is,  the  reliability  of  the  component,  then 

R(t  +  At)  =  R(t)  (1  -  XAt)  +  o(At)  (22) 

which  yields 


R(t  +  At)  -  R^t^  ^  -\R(t)  (23) 

At 


and  in  the  limit,  the  solution  of  the  resulting  differential  equation  is 

R(t)  =  exp(-\t)  (24) 

If  T  is  the  random  variable  denoting  time  to  failure  of  the  component,  then 

ft:  (T  >  t)  =  R(t)  =  exp  (-Xt)  =  l-F(t)  (25) 

where  F(t)  =  l-exp(-Xt)  is  the  exponential  distribution  characterizing  the 
constant  failure  rate,  X. 

If  it  is  assumed  that  the  time-to-failure  distributions  for  individual 
components  are  exponential,  the  best  estimate  of  system  failure-rate  is 
obtained  by  adding  the  component  failure-rates,  and  then 
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RgCt)  =  _TT  R^(t)  =  exp  [  -{X^+X^+   ...+X^)t  J  =  exp(-Xgt)        (26) 

where  R  (t)  is  the  system  reliability 

R. (t)  is  the  reliability  of  the  i-th  component 
"K.  is  the  failure-rate  of  the  i-th  component 

s     is  the  failure-rate  of  the  system 
and    n     is  the  number  of  components  in  the  system 

Epstein  (i960)  enumerates  two  procedures  for  estimation  from  life 
test  data  when  the  underlying  probability  density  function  is  exponential. 
In  the  first  case,  life-testing  is  discontinued  after  a  fixed  number,  say  r, 
of  the  components  have  failed.  In  the  second  case,  the  testing  is  discon- 
tinued after  a  fixed  amount  of  total  life  T  has  elapsed.  In  both  cases,  the 
components  under  test  may  or  may  not  be  replaced  in  case  of  failure. 

Equation  (25)  gives  the  exponential  distribution,  characterizing  the 
constant  failure-rate,  as 

F(t)  =  1  -  exp  (-Xt)  (27) 

The  probability  density  function  for  the  above  distribution  is 

f(t,\)  =  X  exp  (-Xt)  t  >  0,  X  >  0  (28) 

which  in  terms  of  the  mean  life,  e,  is 

f(ti8)  =  -^  exp  (-  |)  t  >  0,  e  >  0  (29) 
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Case  I;  If  the  number  of  components  initially  on  test  is  n  and 

testing  is  discontinued  after  r  items  fail,  then  the  best  estimate,  Q   , 

'  '  r,n' 

of  the  mean  life,  9,  is 


^n=   T^  (30) 


where  T    is  the  accumulated  test  time  until  the  r-th  failure  occurs. 
r,n 

If  the  components  were  not  replaced,  in  case  of  failure,  then 


^r,n=  iil^^("-H,  (31) 


where  t.  (^■t._|_^)  is  the  observed  time  for  the  i-th  failure.  If  the 
components  were  replaced  in  case  of  a  failure,  then 


^r,n  =  "^r  (32) 

The  probability  density  function  of  9   ,  in  either  the  replacement 

r,n  *^ 

or  non-replacement  case,  is 


1    /r\    r-1     /    ,„v 

y    exp  (-ry/8;        y  >  o 


J  (r-D!  (?j 


f^(y)  =V  ^^-^^'-    W  (33) 

elsewhere 


From  quation  (33),  it  follows  that  the  quantity 


2r  9        PT 

^>"  =  -kill  (34) 
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is  distributed  as  chi-square  with  2r  degrees  of  freedom,  and  hence  the 
two-sided  lOO(l-a)  per  cent  confidence  interval  for  9  is 


2T  2T 

^>"  <  Q  <     ,^^'"  (35) 

y^   (2r)  V   (2r) 

7  2 


A  one-sided  lOO(l-Qf)  per  cent  confidence  interval  for  9  is 


e  >  -^ —  (36) 


^  (2r) 


It  is  also  interesting  to  estimate  another  quantity,  t  ,  v/here 

t  is  that  life  such  that 
P 


Pr  (T  >  tp)  =  p  (37) 


For  the  exponential  distribution 


tp  =  ein  (1/p)  (38) 


and  the  maximum  likelihood  estimator  of  t  is 

P 

^p  =  i,n  ^"  ^^/P)  (39) 
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The  lOO(l-a)  per  cent  confidence  intervals  for  t  are 


2T   ln(l/  )  2T   ln(l/„) 

4  (2r)  X?  a  (2r) 

2  2 


(40) 


for  the  two-sided  case,  and 


2T   In(lA) 
tp>    ^i"    /P^  (41) 

4(2r) 


for  the  one-sided  case. 

Case  II;  The  number  of  components  initially  on  test  is  n  and  the 
test  is  discontinued  after  a  fixed  amount  of  total  life  T  has  elapsed.  If 
r  be  the  number  of  items  which  fail  in  the  interval  (o,T),  then,  the  two- 
sided  lOO(l-a)  per  cent  confidence  interval  for  9  is 


2"^  <  e  <      2T      ■  (42) 


2  ^ 

y-a  ^^^2)  Xi  «  (2rf2) 

7  ^"2 


The  one-sided  lOO(l-a)  per  cent  confidence  interval  for  Q  is 


e>  -T^- (43) 

xl  (21^2) 
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The  one-sided  100(1-q?)  per  cent  confidence  interval  for  t  is 


t  >    ^^  ^"  (yp)  (44) 

p  ^  ' 

X^  (2r^-2) 


Zelen  (l96l)  discusses  some  of  the  pitfalls  associated  with  the  use 
of  the  exponential  model  in  life  test.  However,  the  exponential  model  has 
come  to  have  a  special  significance  in  system  reliability  studies  because 
of  the  ease  and  simplicity  with  which  the  prediction  is  achieved. 

The  Weibull  Afodel.  In  the  development  of  the  exponential  model, 
it  was  stated  that  the  failure-rate  of  mechanical  and  electrical  components 
is  distributed  over  time  as  shown  in  Figure  3.  In  practice,  if  a  component 
or  a  system  is  unreliable  in  its  first  hours  of  life,  it  will  be  even  less 
reliable  in  its  hundredth  or  thousandth  hour.  The  failure-rate  curve  of 
Figure  4  indicates  that  the  failure-rate  is  a  monotonically  increasing 
function  with  time  and  that  the  failure-rate  does  not  remain  constant  over 
any  appreciable  duration  of  time. 

If  the  failure-rate  is  a  function  of  the  test  duration  or  the  oper- 
ational age  of  the  component,  then,  equation  (23),  in  the  limit,  becomes 

~-  [R(t)]  =  -\(t)  R  (t)dt  (45) 

the  solution  of  which  is 

R(t)  =  exp  [-H(t)]  (46) 


Time(hours) 


Figure  4:  Failure-rate  versus  test  duration  for  the  Weibull  model 
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where 


H(t)  =  J'  X(r)dT 


(47) 


The  idea  of  a  non-constant  failure-rate  function  is  associated  with 
Weibull  who  showed  that 


H(t)  =  (kt)^ 


(48) 


where  k  and  3  are  positive  real  constants.  Various  values  of  these  parameters 
result  in  distributions  representing  different  types  of  failure  phenomena. 
For  example:  if  3  is  equal  to  1,  the  resulting  probability  of  survival, 
R(t),  exhibits  properties  typical  of  random  failure.  If  3  is  greater  than  1, 
then,  R(t)  exhibits  properties  typical  of  wear-out  failures. 

A  simple  Weibull  cumulative  distribution  function  is  defined  as 


F(t)  =  1  -  exp  \-{t-yf/a   ]      t^Y»a>o,  3>o 


(49) 


where  or  is  the  scale  parameter 

3  is  the  shape  parameter 

Y  is  the  location  parameter 
The  Weibull  probability  density  function  is,  then, 


f(t)  =  Jlkll 


3-1 


a 


exp  1^ 


-(t-Y)%  ] 


(50) 


Electronic  equipment  failures  can  be  classified  as  catastrophic  or 
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sudden  failures  and  wear-out  or  delayed  failures.  For  the  catastrophic 
failures,  a  Weibull  distribution  with  location  parameter  equal  to  zero  and 
shape  parameter  less  than  one  is  an  appropriate  mathematical  model.  In  case 
of  wear-out  failures,  very  few  failures  of  this  nature  occur  before  a  certain 
period  has  elapsed  after  v/hich  the  failure-rate  increases  rapidly  with  time. 
In  view  of  these  considerations,  a  Weibull  distribution  with  some  positive 
location  parameter  and  shape  parameter  greater  than  one  would  be  a  reason- 
able model  Kao  (1959)  has  set  up  another  useful  model  which  is  a  combin- 
ation of  the  above  two  models,  thereby  taking  into  account  both  the  me- 
chanical strength  of  the  item  and  its  wear-out  quality. 
The  Weibull  model  for  catastrophic  failure  is 

F^(t)  =  1  -  exp  (-t^l^,p       t  >  0,  ff^  >  0,  0  <  3i  <  1    (51) 

The  Weibull  model  for  wear-out  failures  is 

F2(t)  =  1  -  exp  [-(t-Y2)'^yof2]   ^  >  Y2»  «2  ^°'  ^2  ^  ^      ^^^^ 
Then,  the  combined  Weibull  model  is 

F(t)  =  pFi(t)  +  qF2(t) 

=  1  -  p  exp  (-t^Va;^)  -  q  exp  [-(■t-Y2)V''2J  C^^) 

where  p  and  q  are  the  proportions  of  catastrophic  and  v^ear-out  failures, 
respectively;  and  p  +  q  =  1, 
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A  graphical  method,  as  presented  by  Kao  (1959),  for  estimating  the 
various  parameters  from  life  tests  data  is  discussed  below: 

The  experimental  cumulative  density  function  of  the  failure  data 
is  plotted  on  the  Vi/eibull  probability  paper*  and  a  VJeibull  plot  is  fitted 
to  this  data,  as  shown  in  Figure  5.  The  tanget  lines,  drawn  at  each  end 

A         A 

of  the  Weibull  plot,  give  the  estimates  pF^  and  qF„  of  pF^(t)  and  qF^(t), 

A  A 

respectively.  The  intersection  of  qF^  with  the  time  scale  gives  Yo  which 

A 

estimates  y_.  From  the  intersection  of  qF„  with  the  upper  borderline  of 

A 

the  graph  paper,  drop  a  vertical  line  whose  intersection  with  pF.  as  read 
on  the  per  cent  failure  scale  gives  p  which  estimates  p  and  hence,  q(=  1-p). 
From  these  estimates  of  p  and  q,  the  combined  Weibull  population  can  be 
separated  into  the  two  sub-populations.  The  estimated  number  of  components 
that  belong  to  F  (t)  and  Fp(t)  are  np  and  nq,  respectively;  and  hence,  the 
cumulative  distribution  functions  F^(t)  and  Fp(t)  can  be  calculated.  The 
y-intercept  and  the  slope  of  the  Weibull  plots  of  the  separated  sub-populations 
F. (t)  yield  the  estimated  of  Ina.  and  3.,  respectively  (i  =  1,2), 

The  combined  Weibull  model  can  be  approximated  by  a  composite  model, 
in  case  of  small  p  and  large  Yo*  The  Weibull  plot,  then,  consists  of  two 
distinct  linear  portions  which  are  treated  as  two  sections  of  a  composite 
population,  as  shown. 


?Jt)   =  1  -  exp  (-t  7^3)        t  ^  6,  ^3  ^i  0,  0  <  ^3  <  1   (54) 


*Weibull  probability  paper  contains  In  versus  In-ln  scales  and  other 
scales  calibrated  for  life-testing  use.  It  is  available  from  Cornell 
University,  Ithaca,  New  York, 
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04/ 
F^Ct)  =  1  -  exp  (-t  1 01^)  t  s  6,  a^  >  0,  3^  >  1      (55) 

where  6  is  the  time  at  which  the  proportion  of  catastrophic  failures  is 
equal  to  that  of  wear-out  failures,  that  is 


1  -  exp  (-6  %^)  =  1  -  exp  {'^''Jot^  (56) 


and  hence 


,.,\^aV'''^-H..,['z^\T\  („, 


3.  -  3' 


The  estimation  of  the  parameters  in  the  composite  model  case  is 
simplified,  as  discussed  below: 

The  experimental  cumulative  density  function  is  plotted  on  a  Weibull 
probability  paper  and  two  straight  lines  are  approximated  among  the  points. 
The  abscissa  of  the  intersection  of  these  two  lines  gives  the  estimate  of  6. 
The  y-intercept  and  slope  of  the  lines  yield  the  estimates  of  Inor,  and  3., 
respectively  (i  =  3,4). 

The  combined  Weibull  model  technique  can  be  extended  to  the  prediction 
of  failure-rate  of  a  complex  system  comprising  of  different  components. 
The  Weibull  cumulative  density  function  for  the  complex  system  is 


n 
F(t)  =   E  p^  F  (t)  0  ^  p.  ^  1  (58) 

i=l 


1 

-1 
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where 

Fj_(t)  =  1  -  exp  [-  (t- 

V.)  ^i/.J 

1 

(59) 

1 

is  the  i 

-th  component  cumulative 

density  function  form  and  p. 

is  the 

1 
pro-           i 

portion 

of  the  i-th  component  in 

the  system 

,  so  that 

n 

2  Pi 
i=l 

=  1. 

\ 

The 

Weibull  probability  den 

sity  function  for  the 

system 

is 

\ 

■ 

f(t)  =  \    p  f  (t) 
i=l  ^  ^ 

i 

(60) 

with  the 

i-th  component  probabil; 

Lty  density 

function 

1 

0.(t-Yi)^i-l 
f.(t)  -  i'         exp  -(t-Yi: 

Vi] 

(61)          ' 

J 

3 

The  reli 

able  life,  9  ,  for  any  specified  r, 

such  that 

i 
i 
1 

r 

(62)          • 

s 

1 

is  given 

by 

e^  =  Y  +  or  "^^  (-Inr)  ^^ 

\ 

1 

1 

:l 

1 

(63)          ! 
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The  Markovian  Model.  The  model  developed  by  Tainiter  (1963)  takes 
into  account  failures  occurring  due  to  drifting  of  the  individual  circuit 
components  from  their  nominal  values  and  also  failures  due  to  the  cata- 
strophic failure  of  any  single  component.  This  model  assumes  that  the  state 
of  a  parameter  of  a  component  at  some  future  time,  given  its  .state  at  the 
present  time,  is  independent  of  the  present  age  of  the  component,  the  number 
of  previous  states  and  the  duration  of  these  previous  states.  It  is  also 
assumed  that  a  change  in  the  value  of  a  parameter  of  one  component  does  not 
change  the  environment  in  which  the  remaining  components  are  operating. 

The  state  of  the  parameter  of  a  component,  at  some  given  time  of  obser- 
vation, is  the  number  assigned  to  the  set  which  contains  the  value  of  the 
parameter.  It  is  assumed  that  the  set  of  states  of  the  system,  which  are 
formed  from  the  component  states,  can  be  partitioned  into  operative  and 
failure  sets  so  that  a  known  functional  relationship  exists  between  the  com- 
ponent values  and  the  output  parameters  of  the  system. 

It  is  assumed  that  a  transition  of  the  system  from  one  state  to  another 
is  merely  a  transition  of  a  single  component  parameter  from  one  state  to 

another  and  that  the  probability  that  two  component  parameters  change  state 

2 
in, a  given  time  interval,  At,  is  the  order  of  (At)  .  From  the  theory  of 

continuous-time  Markov  processes  (Rosenblatt,  1962)  it  can  be  shown  that 

P(t)  =  exp  (\t)  (64) 


where  P(t)  is  a  matrix  whose  elements  P.  .(t)  are  defined  as  the  probability 
that  the  system  is  in  state  i  at  time  zero  jind  state  j  at  time  t;  X  is  a 
matrix  whose  elements  X.  .  are  defined  so  that  X.  .At  +  o(At),  (i?^j),  is  the 
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probability  that  a  system  in  state  i  at  time  zero  enters  state  j  by  time  At 

viith  \. .   =       E  \.  ,,  The  matrices  P  and  X  can  be  obtained  from  the  matrices 

p  and  X  for  the  component  parameters  of  the  system, 
c     c 

Brender  and  Tainiter  (l96l)  have  derived  a  formula  for  computing 
the  reliability  of  a  system;  namely 


R  (t)  =  exp(X^  t).  l"  (65) 

where  R(t)  is  a  vector  whose  elements  R. (t)  are  defined  as  the  probability 
that  the  system  has  not  entered  the  failure  set  of  states  in  the  time  interval 
(o,t)  given  that  the  system  was  in  state  i  at  time  zero.  X  is  obtained  from 
X  by  deleting  those  rows  and  columns  of  X  which  correspond  to  failure  states 
for  the  system.  I  is  a  unit  column  vector  of  conformable  dimensions.  If 
P(o)  is  the  initial  probability  vector  for  the  system,  then,  its  elements 
P.(o)  are  the  probabilities  that  the  system  is  in  state,  i  initially,  i  ranging 
only  over  the  operable  states.  The  inner  product  of  P(o)  with  R(t)  will  yield 
the  unconditional  system  reliability. 

Several  problems  become  apparent  when  it  is  attempted  to  utilize  this 
model  to  predict  the  reliability  of  an  actual  system.  It  is  necessary  to 
define  as  to  what  constitutes  system  failure.  It  must  be  possible  to  obtain 
life-test  data  for  the  time  variations  of  each  of  the  significant  component 
parameters  of  the  system  measured  at  equal  time  intervals.  The  reliability 
prediction  for  the  system  is  based  upon  the  joint  effect  of  the  time  vari- 
ations of  each  component  of  the  system  as  observed  in  life-test  measurements 
of  the  behaviour  of  the  component.  Hence,  the  failure  boundaries  of  the 
circuit  must  be  determined  in  terms  of  the  values  of  the  component  parameter. 
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The  number  of  system  states,  and  therefore,  the  numbers  of  states  for  each 
component  parameter,  must  be  kept  as  small  as  possible  and  it  must  be  de- 
cided with  a  reasonably  small  error  whether  or  not  the  system  is  operable, 
if  the  component  parameters  take  any  value  within  the  ranges  of  the  states. 

The  transition  probabilities  for  each  component  parameter  are  esti- 
mated from  the  definition  of  component  parameter  states.  Let  n^.  denote 
the  number  of  transitions  a  component  parameter  makes  from  state  i  at  a 
given  time  point  to  state  j  at  the  next  time  point.  Let  n.  be  the  total 
number  of  transitions  made  by  a  component  parameter  out  of  state  i.  Then, 

A 

the  maximum  likelihood  estimates  P. .  of  P..  are 


P.  .  =  !ii  (66) 

iJ    n 


The  matrix  X   can  be  obtained  from  the  knowledge  that  a  transition  from 
system  state  i  to  system  state  j  is  represented  by  a  change  of  state  for 
only  a  single  component  parameter.  If  the  change  in  system  from  state  i 
to  j  corresponds  to  a  component  transition  k  to  1,  then 

A 

where  X,   (y)  is  the  k-1  entry  of  the  matrix  X^  (y)  for  the  component  y. 
Thus,  the  system  matrix  X  will  have  a  large  predetermined  number  of  zeros, 

A 

which  also  simplifies  the  reliability  calculation.  The  estimate,  X  ,  of  the 

A 

component  matrix  X  is  obtained  from  the  estimated  matrix  P  for  the  component 
c  c 
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from  the  follovdng  equation 


\^   =  log  P^  (68) 


A 

Equation  (68)  is  useful  only  when  log  P  converges;  for  which  P. .  must 

c  ij 

be  greater  than  0.5,  for  all  i.  This  condition  seems  to  be  adequate, 
because  for  a  reliable  design,  the  probability  that  the  component  does 

A  A 

not  change  state  is  greater  than  one-half.  Having  obtained  X  and  thus  X, 

c 

the  system  reliability  iscomputed  with  the  use  of  equation  (65)» 

The  Markovian  model  provides  an  analytic  method  of  computing  system 
reliability  from  component  drift  and  failure  data.  The  other  models  for 
computing  system  reliability  define  reliability  as  the  probability  of 
failure  at  a  given  time.  The  Markovian  model  goes  a  step  further  in  con- 
sidering reliability  as  the  probability  of  failure  in  an  interval  of  time. 
The  model  is  not  useful  for  highly  complex  systems  for  which  the  failure 
characteristics  are  described  by  equations  whose  solutions  require  an 

exorbitant  amount  of  time  even  on  a  high  speed  computer. 

The  VJorst-Case  Design  Model  (Combs;  1963).  A  system  is  said  to 

be  worst-case  designed  if  and  only  if  it  performs  adequately  vdth  the 
component  parameters  assuming  values  at  every  combination  of  their  combi- 
nation limits.  If  the  output,  y,  of  a  system  be  a  function  of  its  input 
and  its  component  parameters,  then 

y  =  f  (xj^,  x^,    ...,  x^)  (69) 


where  x.  represents  either  an  input  or  a  component  parameter.  If  (x.) 


max 
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and  (x^)^^j^  are  the  maximum  and  minimum  tolerance  limits,  respectively, 
for  the  component  x^,  the  2"  possible  values  of  y  from  equation  (69) 
also  have  a  maximum  and  a  minimum  value.  If 

y    ^  maximum  acceptable  output 


'max 


and 

^min  ^  "iininiuni  acceptable  output 

then,  the  system  is  worst-case  designed. 

Let  X  be  the  failure-rate  of  a  component  having  a  cumulative  failure 
distribution  function  of  the  form  1-  exp  (-Xt).  If  p  and  (l-p)  be  the 
proportions  of  drift  failure  and  catastrophic  failures,  respectively,  then 

pX  =  drift  failure-rate 
(l-p)X  =  catastrophic  failure-rate  (7O) 


If 


Fq  (t)  =  1  -  exp  (-pXt) 

(71) 
and  F^  (t)  =  1  -  exp  [-(l-p)Xt] 

be  the  drift  and  catastrophic  failure  distributions,  respectively,  and  if 
the  drift  and  catastrophic  failures  be  assumed  to  be  independent  events, 
then  the  component  failure  distribution  is 
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F(t)  =  F_(t)  +  F  ft)  -  F^(t)F  (t) 


D 


D^"'  c' 


=  [l-e-P^tJ  ^  ^^_^-(i.p)xtJ  _  ri_3-pU-]  ri.,-(l-p)Xt-]    (72) 


For  a  system  with  n  components,  the  failure  distribution  of  the  system 


is 


n     i=i  I-       J   L   1=1  J 


y.  (-e-^^^^J  [:-J^  e-'-Pi)Ht] 


(73) 


Consider  a  component  with  the  highest  failure-rate  X  for  this  system. 
If  the  system  is  assumed  as  being  made  up  of  n  of  these  components,  then, 
the  failure  distribution  for  a  worst-case  designed  system  is 


F^(t)  =  e 


-n(l-p)X^t  /   -PX  tvn    -n(l-p)\  t 


(l-e-^^™^) 


-••  e 


m 


+  1 


(74) 


=  1  + 


,  k  m 


(75) 


(^, 


where  c^  =  kp  +  n(l-p)  and  {  !}   are  binomial  coefficients 
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The  probability  density  function  for  the  system  failure  is 


h]  (-1)^'"''''"*  (-vj  (^*) 


The  mean-time-to-failure  for  the  system  is 


e=rtf(t)dt=  E  i-D^-'U'}]    -i-  (77) 

N  /     km 


For  the  special  case  of  p=o,  that  is,  no  drift  failures 


e  =  _1_  (78) 


nX 
m 


and  for  perl,  that  is,  no  catastrophic  failures 


9=    1    2   i  (79) 

\       k=l  k 


Thus,  the  simplified  case  of  the  worst-case  design  model  yields 
equations  (78)  and  (79)  which  appear  frequently  in  the  reliability  liter- 
ature as  the  equations  of  mean-time-to-failure  of  a  series  and  parallel- 
redundant  system  of  n  components,  respectively. 
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CORRELATION  TECHNIQUES 

In  this  section,  two  analytical  treatments  are  discussed  which  indi- 
cate the  existence  of  a  definite  correlation  between  the  raean-time-to-failure 
under  normal  operating  conditions  and  that  under  accelerated  environmental 
conditions.  The  results  of  these  two  treatments  prove  the  validity  of  our 
assumptions  that  true  acceleration  must  be  achieved  and  that  the  acceleration 
factor  must  not  be  too  large  so  as  to  introduce  a  change  in  the  mode  of 
failure. 

The  analytical  treatment  is  restricted  to  the  combination  of  only  two 
types  of  environments,  E  and  E^,  say.  Let  E.  .  and  E„  .  represent  the 
i-th  and  j-th  stress  levels  of  the  environments  E  and  E^,  respectively. 
Let  e.  .  denote  the  mean-time-to-failure  of  a  system  under  a  particular  combi- 
nation  of  the  i-th  stress  level  of  E.  and  the  j-th  stress  level  of  E^. 

If  it  is  assumed  that  the  mean-time-to-failure  is  a  negative  exponential 
function  of  the  environmental  stress  level,  then 

e.  .  =  a  exp(-b  E   )  i  =  1,  2,  ...  (80) 

for  a  particular  value  of  E„  .  (j  =  1,  2,  ...).  The  constants  a  and  b.  are 
the  intercept  on  the  mean-time-to-failure  axis  and  the  slope,  respectively, 
of  the  semi-logarithmic  plot  of  the  mean-time-to-failure  curve s^  A  set  of 
these  curves  can  be  obtained  by  varying  the  stress  level  of  the  environment 
E  over  its  entire  range,  for  a  particular  level  of  £„;  and  then  repeating 
this  for  all  other  values  of  Ej. 

For  two  different  values  of  j,  as  shown  in  Figure  6,  the  correlation 
coefficient  between  the  two  mean-time-to-failure  curves  is 
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fa,,.,=  if>Ma/¥l.i,[3..e-V<^i^^^>laH,.,    (31) 


^l.i 


where 


El,  i' 

6E  =  f  '   dE.  .  =  E.  .  ,  -E,  .  (82) 

Jc      1,1    1,1    1>  1 

1»  1 


is  the  separation  parameter  which  is  related  to  the  acceleration  factor. 
Equation  (81)  yields 


e(J,  j',  6E)  =  K(j,  j')  e"^j'^^  (83) 

6E 


The  graph  of  normalized  correlation  coefficient  versus  the  separation 
parameter,  see  Figure  7,  indicates  that  for  reasonable  correlation  to  exist, 
the  separation  parameter  and  therefore,  the  acceleration  factor  should  not 
be  too  large. 

If  the  mean-time- to- failure  is  not  a  simple  negative  exponential  function 
of  the  environmental  stress  level,  as  shown  in  Figure  8,  it  can  be  divided 
into  different  sections,  such  that,  each  section  has  its  own  characteristic 
exponential  behaviour,  that  is 
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Figure  6:  Mean-time-to-failure  versus  environment 
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Figure  7:  Auto-correlation  function  versus  lag 


Figure  8:  Semi-logarithmic  plot  of  mean-time-to-failure 
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a.  exp  (-b.E,  .) 


e.  .  =  <a' .  exp  (-b'  .E,  .) 


a",  exp  (-b'.E,  . ) 


o  ^  i  ^  i' 


i'  ^  i  S  i"  (84) 


.  ^  .11 


for  a  particular  E„  .,  However,  in  order  to  achieve  any  useful  prediction 
of  the  system  mean  life  under  normal  environmental  conditions  from  acceler- 
ated life-test  data,  it  is  essential  that  the  points  of  discontinuities  of 
the  slope  of  the  curve  be  known.  If  the  accelerated  life-tests  yield  data 
in  the  region  beyond  the  stress  level,  E.  .„,  then,  the  extrapolation  of 
the  curve  assuming  the  same  exponential  variation  with  environment  would 
yield  an  error  in  the  prediction  of  mean-time-to-failure  under  normal  en- 
vironmental conditions,  as  shown  in  Figure  8, 

The  anomalies  brought  out  in  the  above  discussion  could  be  irradicated 
if  the  variation  of  the  failure-rate  of  the  system  with  environmental  stress 
level  is  assumed  to  be  of  the  form 

\^^.  =  I   -  exp|_-(E^^^)  ^/^.   J       i  =  1,  2,  3,...         (85) 

for  a  particular  E„  .  (j  =  1,  2,  ,,,),  Equation  (85)  is  similar  to  the 
Weibull  probability  density  function  and  the  parameters  a.   and  15.  can  be 
estimated  in  the  same  way  as  discussed  in  the  section  on  the  Weibull  model. 

This  technique  yields  good  results  in  cases  where  the  acceleration 
curve  is  smooth  and  well  defined  over  a  range  of  environmental  stress  levels 
and  no  sharp  discontinuities  appear  at  or  around  the  stress  level  at  which 
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the  failure-rate  is  to  be  predicted. 

It  has  been  discussed  in  the  sections  on  the  life-testing  techniques 
and  the  exponential  model  that  the  estimation  of  system  failure-rate  is 
subject  to  variation  within  the  confidence  interval  (dotted  lines),  in 
Figure  9.  The  width  of  these  confidence  interval  depends  on  the  degree 
of  confidence  placed  in  the  estimation,  A  l-or  per  cent  confidence  inter- 
val, CI,   ,  implies  that  the  true  percentage  of  the  observations  lie  within 
the  confidence  interval  states,  subject  to  a   per  cent  error. 

The  system  failure-rate  for  a  set  of  environmental  conditions  is,  say, 
normally  distributed  about  its  mean  value.  The  l-or  per  cent  confidence 
interval  is  then  obtained  by  the  use  of  statistical  procedures.  The  confi- 
dence interval  for  the  associated  correlation  of  the  data  under  two  different 
sets  of  environmental  conditions  may  also  be  similarly  calculated.  The  latter 
part  of  this  section  is  devoted  to  the  development  of  the  correlation  coef- 
ficient between  two  random  variables,  namely,  the  failure-rates  under  dif- 
ferent sets  of  environmental  conditions. 

The  degree  of  variation  in  the  failure-rate  of  the  system  corresponding 
to  a  given  stress  level  of  environment  E.,  measures  the  closeness  of  the  corre- 
lation. In  the  present  discussion,  it  is  assumed  that  for  the  correlation 
between  the  system  failure-rate  and  the  stress  levels  of  the  environment  E  , 
E-  can  be  maintained  constant  at  its  j-th  stress  level.  Usually  there  is  no 
means  of  varying  the  environment  E  quickly,  v;hile  E^  remains  constant,  in 
order  to  observe  the  effects  upon  the  system  failure-rate.  In  the  latter 
case,  the  problem  is  one  of  multiple  or  partial  correlation.  However,  the 
present  discussion  deals  only  with  the  correlation  of  two  variables,  >.  and 
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E^,     Let 


\    =  a.   exp(b.  E.  .)  i  =  1,  2,  ...         (86) 

s    J      J  J-  >i 


for  a  particular  E^  .  (j  =  1,  2,  ),  Hence 


in  X^  =  m  a.  +  b.  E^^.  (87) 


Replacing  In  X  by  Y  and  E  .  by  X 


Y  =  In  a.  +  b.X  (88) 


Using  Pearson's  formula  (Chaddock,  1925)  the  coefficient  of  correlation  is 


p  =  ^y  (89) 

''   ^Vy 


where  x  and  y  are  the  deviations  of  each  X  and  Y  from  their  respective  means, 
a„  and  a^   are  the  standard  deviations  of  the  entire  X  and  Y  distributions, 
respectively,  and  N  is  the  total  number  of  related  pairs  of  observations. 
If  a  very  large  number  of  similar  random  samples  of  the  same  system 
are  examined  and  related,  P  is  shown  to  be  itself  a  variable.  The  signifi- 
cance of  ^  evidently  depends  upon  the  amount  of  this  probable  variation 
which  is  ue  to  the  uncontrolled  conditions  of  sampling. 


Pr  (error  in  ^ )  =  0,6745  LjlSL  (90) 
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Conservative  statistical  practice  in  interpreting  ^  requires  that  the  size 
of  the  coefficient  should  be  at  least  four  times  its  probability  of  error 
before  it  becomes  indicative  of  any  significant  degree  of  correlation.  It 
has  been  strongly  emphasized  that  a  ^  of  low  value  to  be  significant  must 
be  based  upon  many  more  cases  than  one  of  high  value. 

The  coefficient  of  correlation  is  a  pure  number  indexing  the  degree  of 
correlation  betv;een  two  variables.  It  does  not  indicate  how  much  variation 
on  the  average  may  be  expected  in  one  variable  when  the  amount  of  variation 
in  another  related  variable  is  known.  However,  if  the  regression  equation 


a. 


Y-Y=^-I  (X-X)  (91) 


^X 


is  considered,  it  gives  a  more  complete  description  of  the  relationship  than 
the  coefficient  of  correlation. 

The  predicted  single  value  of  Y,  from  equation  (91 ),  corresponding 
to  a  given  value  of  X  is  subject  to  a  range  of  error  determined  by  the 
scatter  about  the  straight  line  of  average  relationship,  as  shown  in  Figure  10. 
The  less  this  scatter,  the  more  closely  it  approaches  the  condition  where 
there  would  be  only  one  value  of  Y  for  a  given  value  of  X.  The  amount  of 
scatter,  S^.,  in  the  Y  variable  is  given  by 


^Y  =  ^Y  \J'-V 


(92) 


If  the  distribution  of  Y  is  assumed  to  be  normal,  then,  in  about  68  per  cent 
of  the  cases,  the  actual  values  of  Y  will  not  differ  more  than  ±  S  from  the 
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X 


Figure  10:  Distribution  of  the  error  about  the  regression  line  RR 
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predicted  values, 

Pearson's  product-deviation  method  of  measuring  the  coefficient  of 
correlation  P  between  two  variables  is  based  upon  the  hypothesis  that  a 
straight  line  fits  most  closely  the  means  of  Y  and  X  in  a  correlation 
table.  But  sometimes,  the  means  conform  more  closely  to  some  other  form 
of  curve.  When  the  line  of  the  means  is  non-linear,  a  low  of  value  of  ^, 
as  calculated  from  the  straight  line  fit,  does  not  necessarily  indicate 
that  the  degree  of  correlation  is  really  small  or  that  the  two  variables 
are  unrelated. 

For  measuring  non-linear  relationships  Pearson  has  proposed  a  new 
constant,  Tj  ,  the  correlation  ratio. 


11   -  cr  of  the  means  of  Y  for  a  particular  X 
'yX~  a  ot  the  entire  Y  distribution  ot  ail  X 


that  is 


Z:n^(Y.,  -  Y)^ 


\x  = 5^ (^^) 


where  n^  is  the  total  frequencies  for  a  single  value  of  X 

Yy.  is  the  mean  of  Y  corresponding  to  a  particular  value  of  X 

Y  is  the  mean  of  the  entire  Y  distribution  of  all  X 
N  is  the  total  number  of  related  items  in  table, 
and   o"y  is  the  standard  deviation  of  the  entire  Y  distribution. 
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The  correlation  ratio  is  a  measure  of  the  degree  of  correlation 
applicable  to  both  linear  and  non-linear  relationships.  It  is  limited 
in  that  it  is  not  applicable  to  ungrouped  data  and  in  that  it  does  not 
enable  one  to  estimate  values  of  the  related  variable  from  known  values 
of  the  given  variable,  as  in  the  case  of  regression  equations  of  the 
straight  lines.  Its  value  lies  in  giving  an  index  of  maximum  correlation. 
It  furnishes  a  means  of  detecting  divergence  from  linear  relationship, 
and  its  use  prevents  errors  in  conclusions  due  to  the  wrong  assumption, 

RELIABILITY  PREDICTIONS  FOR  A  SIMPLE  AMPLIFIER 

Consider  a  simple  amplifier  circuit,  as  shown  in  Figure  11,  The 
failure-rates  of  the  resistors,  capacitors  and  tubes  over  a  range  of  ambi- 
ent temperature  and  for  a  particular  stress  ratio  of  applied  to  rated  wattage 
are  obtained  from  Military  Standardization  Handbook  (MIL-HDBK-217),  With 
the  help  of  techniques  discussed  in  the  handbook,  the  component  failure- 
rates  are  calculated  and  listed  as  X^,  X    and  A^  for  the  resistor,  capacitor 
and  tube,  respectively,  in  Table  I. 

Assuming  that  the  failure-rate  of  each  omponent  remains  constant  over 
the  operating  life  of  the  amplifier,  the  system  (complete  amplifier  circuit) 
failure-rate,  using  the  exponential  model  is  obtained  as  the  sum  of  the  com- 
ponent failure-rates,  that  is 

\  =  ^\^^\'-\  (94) 

neglecting  the  failure-rates  of  the  power  source  and  connector  leads,  A  semi- 
logarithmic  plot  of  the  system  failure-rate  versus  temperature  does  not  yield 
a  straight  line  over  the  entire  range  of  temperature.  Hence,  any  effort  to 
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Figure  11:  A  simple  amplifier  circuit 
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Figure  12:  Tube  failure-rate  adjustment  factor  for  heater  voltage 
(MIL-HDBK-217) 
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Figure  13:  Adjustment  factor  for  temperature/dissipation  effects 
(MIL-HDBK-217) 
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Figure  14:  Adjustment  factor  for  temperature/dissipation  effects 
(MIL-HDBK-217) 
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Figure  15:  Failure-rates  for  1V1IL-R-26C  quality  power  wirewound  resistors, 
temperature  range  V  and  Y  (relative  humidity  less  than  60^) 
(MIL-HDBK-217) 
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Figure  16:  Failure-rates  for  MIL-C-5B  quality  mica  capacitor  (molded 

type),  temperature  range  P  (relative  humidity  less  than  60^) 
(MIL-HDBK-217) 
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Table  1.  System  failure-rate  for  the  amplifier  circuit. 


Temperature 

ratio  of  a 

pplied  to 

rated  wattage  =0.5 

^^ 

X 
c 

h 

\ 

10 

0.038 

0.041 

0.130 

0.326 

20 

0.041 

0.044 

0.130 

0.341 

30 

0.045 

0.048 

0.130 

0.361 

40 

0.048 

0.053 

0.130 

0.380 

50 

0.053 

0.058 

0.130 

0.405 

60 

0.057 

0.064 

0.130 

0.429 

70 

0.063 

0.071 

0.130 

0.461 

80 

0.068 

0.079 

0.130 

0.492 

90 

0.075 

0.088 

0.130 

0.539 

100 

0.082 

0.099 

0.130 

0.574 

110 

0.090 

0.112 

0.130 

0.624 

120 

0.099 

0.126 

0.130 

0.679 

130 

0.108 

0.145 

0.132 

0.746 

140 

0.119 

0.165 

0.134 

0.821 

150 

0.131 

0.190 

0.136 

0.909 
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predict  the  system  failure-rate  at  normal  temperature,  of  say  30  C,  from 
a  range  failure-rate  data  over  a  range  of  accelerated  temperatures,  say 
100°  C  to  150°  C  by  means  of  this  semi-logarithmic  plot  would  involve 
considerable  error,  similar  to  that  indicated  in  Figure  8. 

However,  if  the  system  failure-rate  is  plotted  versus  temperature 
on  a  Weibull  probability  paper  it  yields  two  distinct  straight  lines  over 
the  entire  temperature  range.  The  application  of  Berrettoni's  (1964)  tech- 
nique for  evaluating  the  constants  a',  and  13.,  using  the  characteristics  of 
the  Weibull  paper  yields 


X.  .  -  ^ 


0.15/  ^  ,^.  c    <-  Q-^o 


fl  -  exp  (-E^  .^•■'7  3.67)       E    ^  83"  C 
1  -  exp  (-E^^i^'^y  6x10^^)      E^^^  ^  83°  C 


(95) 


for  a  particular  E^  .  (=  0.5).   (Refer  to  Figure  17). 

Using  equation  (95)  to  predict  the  system  failure-rate  at  a  temperature 
of  30°  C,  a  value  of  0.364  failures  per  1000  hours  is  obtained,  which  agrees 
satisfactorily  with  the  observed  system  failure-rate  of  0.361  failures  per 
1000  hours  from  Table  I, 

CONCLUSIONS 

The  rapid  developments  in  military  and  space  systems  of  today,  and  the 
stringent  need  for  adhering  to  specifications  provide  only  one  feasible  way 
of  predicting  the  reliability  of  these  systems  and  that  is,  through  acceler- 
ated life-testing.  The  testing  of  a  component  under  severe  environmental 
conditions  accelerates  its  failure-rate  and  enables  prediction  in  a  short 
period  of  time  from  a  small  statistical  population,  thereby  enabling  con- 
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siderable  saving  in  time  and  cost. 

The  accelerated  failure-rate  curves  for  the  system  are  obtained  from 
failyre-rate  curves  of  its  components  by  means  of  various  reliability 
models.  Among  the  four  models  discussed,  the  exponential  and  the  Weibull 
models  afford  simplicity  in  the  evaluation  of  the  system  failure-rate. 
The  Markovian  model  yields  very  satisfactory  results  provided  the  proba- 
bility matrix  of  component  parameter  drift  from  one  state  to  another  is 
known.  The  worst-case  design  model  is  a  mere  extension  of  the  exponential 
model. 

The  correlation  techniques  which  are  discussed  enable  the  evaluation 
of  system  failure-rate  under  normal  accelerated  environmental  conditions. 
These  techniques  yield  satisfactory  correlation  provided  the  accelerated 
environmental  conditions  are  within  the  limits  of  system  specifications, 
and  provided  the  failure  modes  of  the  system  remain  same  over  all  sets  of 
environmental  combinations. 

The  limited  data  available  in  military  specification  handbook  have  been 
used  to  evaluate  the  failure-rate  of  an  amplifier  circuit.  The  failure-rate 
under  normal  environmental  conditions  as  predicted  from  the  failure-rate 
data  under  accelerated  environment  are  seen  to  be  in  close  agreement  with 
the  observed  data,  . 
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APPENDIX  I 

Let  9  be  the  mean  life  of  the  component,  in  hours; 

T  be  the  mean  life  estimated  from  r  failures  in  a  sample 
of  n  tested; 
and     X  be  the  failure-rate,  in  per  cent  per  1000  hours. 

The  required  reliability  is  99.99  per  cent  per  1000  hours,  in  other 
words,  the  failure-rate  is  0.0001  per  1000  hours,  so  that 

X  =  0.0001   lO"''  failures  per  hour 
1000 

and  the  mean  life  is 

e  =  i  =10*^  hours 

If  the  confidence  coefficient  is 

1  -  a  =  0.99 
that  is 

a  =  0.01 
then,  for  a  tv;o-tailed  chi square  prediction  and  r  =  1 


2r 

2 

_ 

2 

a 

2 

X 

o.ol 

9.21 

2 

2 

and,  hence,  for  99  per  cent  confidence 


T  =  lo"^  X  9.21  w  4.6  X  lo'''  hours 
2 
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ABSTRACT 

Accelerated  life-testing  provides  the  only  economically  feasible  way 
of  predicting  the  reliability  of  a  complex  missile  control  system  compris- 
ing of  thousands  of  components,  wherein  each  component  must  be  highly 
reliable.  If  the  accelerated  failure-rate  curves  of  the  individual  com- 
ponents are  known,  then,  the  exponential  reliability  model  can  be  used  to 
obtain  the  system  failure-rate  as 


n 
X  =  S  X.  (1) 

^   i=l  ^ 


where  X.  is  the  failure-rate  of  the  i-th  component  wdiich  is  assumed  to  be 
constant  over  the  useful  operating  life  of  the  system.  The  reliability  of 
the  complex  system  is,  then 


Rg  (t)  =  exp  (-Xgt)   -  (2) 


For  a  system  subject  to  various  combined  stress  levels  of  two  en- 
vironments, the  correlation  between  failure-rate  data  under  accelerated 
and  under  normal  operating  environmental  conditions  is  obtained  by  assuming 
that 

e^^j  =  a^.  exp  (-bjE^^^)        i  =  1,2,...,  k^  (3) 


for  a  particular  E„  .  (j  =  1,2, . ., ,  k^).  6,  .  is  the  mean-time-to- 
failure  (reciprocal  of  the  failure-rate)  of  the  system  under  a  combi- 
nation of  the  i-th  stress  level  of  environment  E  and  the  j-th  stress 
level  of  environment  E^, 

This  correlation  technique  assumes  that  a.  and  b.  are  constant 

J     J 

for  the  entire  range  of  the  stress  levels  of  the  environment,  E. ;  or  in 
other  words,  the  failure  modes  of  the  system  are  the  same  under  combined 
environmental  stress  levels  (E  . ,  E_  .)  for  a  given  j  and  any  i  within 
the  specified  range  of  E..  In  general,  the  quantities  a. and  b.  are 
functions  of  the  combination  environmental  stress  levels  and  represent 
the  intercept  on  the  mean-time-to-failure  axis  and  the  slope,  respectively, 
of  the  mean-time-to-failure  versus  environmental  stress  level  when  plotted 
on  a  semilog  paper.  If  the  points  of  discontinuities  of  the  slope  of  such 
a  plot  are  not  known,  a  considerable  error  will  be  introduced  in  the  pre- 
diction of  system  failure-rate  under  normal  operating  conditions  from 
accelerated  life-test  data, 

A  more  definite  prediction  can  be  achieved  if  it  is  assumed  that 


X^^j  =  1  -  exp{  -(E^^^)  y^j  ]•         i  =  l,2,...,k^    (4) 


for  a  particular  E„   (j  =  1,2,...,  k^),  where  or.  and  3 ,  are  the  scale  and 
•^>J  •^         3  3 

shape  parameters,  respectively,  and  can  be  obtained  graphically.  The  failure- 
rate  under  normal  operation  conditions  can,  then,  be  evaluated  from  equation 

(4). 
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